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Abstract 

The isotropic limit of spin systems with orbital degeneracy is shown to have 
global SU{A) symmetry. On many 2D lattices, the ground state does not 
posses long range order, which may explain the observed spin liquid proper- 
ties of LiNi02- In the SU{4) Neel ordered state, spin-spin correlations can 
be antiferromagneitc between two neighboring sites with parallel magnetic 
moments. 
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In many transitional metal oxides, the electron configuration on the metal ions has or- 
bital degeneracy in addition to spin degeneracy. In these systems, the sign and magnitude 
of the spin-spin couplings depends on the orbital occupancy. This may result in interesting 
magnetic properties of the Mott insulating phase, and is believed to be relevant to unusual 
properties of many vanadium, titanium, manganese and nickel oxides JT|-|r2||. The Hamil- 
tonian describing spin s = 1/2 systems with a two-fold orbital degeneracy (isospin r = 1/2) 
was derived by Castellani et al. ||. The Hamiltonian is rotationally symmetric in s-space, 
but not in r-space. The aniostropy of the latter is due to Hund's rule and the anisotropy in 
orbital wavef unctions. 

In this Letter we study a simplified Hamiltonian, Eq. (1), where we focus on the rota- 
tionally symmetric part in r-space. The insight learned from this higher symmetric model 
should shed light on our understanding of more realistic systems. 

We show the isotropic Hamiltonian (1) has a S77(4) symmetry, and the spin and orbital 
states are described by four flavors. The simplest SU(4) singlet (flavorless) is a four-site 
cluster, in analogy to the 377(2) singlet of a two-site pair. The model provides a new pos- 
sibility for spin liquid ground states in higher dimensions. For the square lattice, using 
the fermion mean field theory, we find the flavor liquid state to be stable against flavor or 
generalized spin density wave formation. By comparing the energies of long-ranged ordered 
states to short-ranged ones on the triangular lattice, we argue the ground state is likely to 
be a resonant plaquette flavor liquid. In the SU(4) Neel ordered state, the spin-spin corre- 
lations can be antiferromagnetic (AF) between two neighboring sites with parallel magnetic 
moments. 

The simplest quantum spin- 1/2 system with two-fold degenerate orbitals (r-1/2) and 
rotational invariance in both s and f -spaces is given by [[J 

H = Jj2( 2g i ■ *j + 1 /2)(2r 4 • tj + 1/2), (1) 

where (ij) is the nearest neighbor (n.n.) pairs. For the AF coupling, which we consider 
exclusively in this paper, J > 0, and we set J = 1 hereafter. Apparently, ([!]) has SU(2) x 
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SU{2) symmetry, representing rotational invariance in both spin and orbital spaces, and also 
interchange symmetry between spins and orbitals. We now show that the full symmetry of 
(|I|) is actually the higher symmery group SU(A), which unifies the spin and orbital degrees 
of freedom. 

Let us start with an intuitive way to see the SU (4) symmetry by rewriting ([!]) as, 

tf = (V4)£(E^'4 + 1 )> (2) 

where A 1 = 2s a , 2r a , 4s a r 13 for a, (3 — x, y, z. A 1 can be considered as the fifteen generators 
of the SU(4) group. We now write them in terms of the more standard generators of group 
theory. 

The Hamiltonian (|I|) acts on a Hilbert space of 4 basis states at each site. Choosing 
these as | s z , r z ), we label them as 

| 1) =|l/2,l/2), | 2) =|-l/2, 1/2), 

| 3) =|l/2,-l/2), | 4) =|-l/2, -1/2). (3) 

These basis states form a fundamental representation of SU(4). The conventional SU(4) 
generators acts on a basis state | /z) according to \ jj) = 5 Utfl \ m). The S^'s are not 
hermitian, instead (S^y = S™. Among the fifteen generators S^, twelve are non-diagonal 
with six raising and six lowering operators. There is an identity Y^tn=i = 1, so that there 
are three independent diagonal ones. The S^s are related to s and f by, 

n z \ / qm om+l\ + \ nm+1 

z * ~~ Z-^i \°m °m+l)i b ~ °m ) 

m=l,3 m=l,3 

2r z = £ i^-S^i), r+ = £ S™ + \ (4) 

m=l,2 m=l,2 

where = s x ± is y , and r ± = t x ± ir v . The commutation relations for spin and orbital 
operators follow from the SU{^) ones, [S^, S l k ] = S n ,kS l m — 5 m ^S k . In terms of S^, ([!]) 
becomes 

# =£S£(W0'). (5) 

(hi) 
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The repeated indices n, m are summed in Eq. fl5|) and hereafter. We conclude that H has 
global SU(A) invariance. 

Eq. (|I]), or equivalently (|5|), gives the effective Hamiltonian for the corresponding Hub- 
bard model in the large U-limit and at 1/4 filling, 

H H = -t ( C L> c i,M + h - c -) + u H n i}IM ni }iy , (6) 

where Ci iM and c|„ are annihilation and creation operators of an electron at site i and state 
| /i). In terms of electron operators, S^(i) = c' im Ci jn . Eq. (§) is also equivalent to one of a 
class of models that has been solved by Sutherland in 1-dimension (ID) fl3| . 

We remark that the AF SU(4) model here is different from the SU(N) model studied 
by Sachdev and Read [f[4| , and by Arovas and Auerbach [15]. These authors considered AF 



SU (N) model on bipartite lattices, where the two sublattices have conjugate representations 
with respect to each other ("quarks" and " antiquarks" ) . In the present model, all the sites 
have the same representation, which is not self-conjugate. 

To get insight on the physical properties, we first consider systems with a few sites. Since 
H has global SU(4) invariance, the eigenstates are given by irreducible representations of 
SU(4). In Fig. 1, we show the Young tableaux for two- and four-site systems. In the two-site 
system, the lower energy (e = —1) states are 6- fold degenerate (total spin s = 1 and total 
orbital r = or s = and r = 1), and higher energy (e = 0) states are 10-fold degenerate 
(s = t = 1 or both = 0). In the four-site system, the ground state is a unique SU(4) singlet 
| SGL), which is rotationally invariant under the SU(4) generators, 

Y,A]\SGL) = 0, (7) 

i 

In terms of S™, the singlet satisfies £i(S£(») - <W 4 ) I SGL) = 0. A SU{4) singlet is 
a singlet of spin, orbital, and the orbital-spin crossing operator U a,/3 = 4s a r /3 , and is a 
generalization of the SU(2) singlet of spin only systems. The energy of the SU(4:) singlet 
of the four-site is found to be — Nj,, with Nb the number of pairs (ij) in (JI|). Hence, Nb = 4 
for a 4-site ring, Nb = 3 for an open chain, and Nb = 6 for a tetrahedron. It is interesting 



to note that the energy of each bond in the four-site system is e& = — 1, the best energy a 
single bond can have. This would be difficult to understand from the conventional valence 
bond picture for spin only systems, and again indicates the difference between ([]]) and spin 



only models including the four-site plaquette RVB state recently discussed in literature |16 
In terms of the fermion operators of (|5]), the SU(4:) singlet can be written as, 

I SGL) = JL £ c J lC ] 2C t 3C t 4 I o), (8) 



{ijki} 

where the sum is over all the permutations of the four sites ijkl = 1234. 

For large number of sites, N s , the system can be a SU(4) singlet only if N s = AN, 
with N integers. This can be proved easily from (|7|). A necessary condition for a 577(4) 
singlet is J2i sf = J2i T l = J2i sfif = 0. This requires 1/4 of the sites of the system in each 
of the four flavor states | /i). Systems with N s ^ 4N may be considered as "edge state" 
or excitations of AN ones similar to odd-number-site systems in the AF spin-1/2 systems. 
For an AF SU(2) spin-1/2 chain, there is a theorem [17] that the ground state is a unique 



spin singlet for finite N s and, in the thermodynamic limit, is either gap less or has broken 
translational symmetry (dimerization). Affleck and Lieb generalized this to all 1/2-integer 
S, and to SU(N) |1£|. The SU(N) representations for which their proof applies includes the 
fundamental representation, which is the representaion for our problem. Thus, the ground 
state of Eq. (JJ) for a ID chain is a unique SU (4) singlet for finite N s = AN, and in the 
thermodynamic limit is either gapless or breaks translational invariance ("quadmerization"). 
Their theorem for S77(4) can be extended to 2D in the same way as for the SU(2) case, but 
requires a long narrow strip as discussed by Affleck |T9| . 



The SU(4) symmetry identified for model (|TJ) has interesting consequencies. Provided 
there is no symmetry breaking, it follows from the symmetry that the thermodynamic cor- 
relation functions, denoted by (...), 

W> = W> = (Krf %rf) = (9) 
where is a function of i and j, and is independent of the indices a and (3 =x, y, z. This 



symmetry has been observed in quantum Monte Carlo calculations of the ID system |20 



For translational invariant systems, the n.n. correlation is related to the energy per bond, 

w = ^ e »-\)> ( 10 ) 



For a ID chain, Sutherland's Bethe ansatz solution for the equivalent model [T3| gives 



three gapless modes. This may be interpreted in our model as pure spin, pure orbital, and 
orbital-spin crossing excitations. The correlation function Wij is expected to decay in a 
power law at T=0, and its sign to be periodic with positive sign if j — i — AN, and negative 
otherwise. The reason for the latter is the tendency for every four neighboring sites to form 
a SU(A) singlet. 

The nature of the ground state of (|]) is of great potential interest. There have been 
numerous theoretical activities since the discovery of high temperature superconductivity to 
find possible spin liquid ground states in two or three dimensions. The additional orbital 
degrees of freedom provides a new possibility for such states. Since there are four equivalent 
single site states | //) in ([I]) in comparison with two states in the spin only systems, we 
expect quantum fluctuations to be stronger, making it more difficult to establish long-range 
order, hence favoring flavor liquid states. 

To illustrate this, we consider model ([!]) on a square lattice and carry out a fermion 
mean field theory. In fermion representation, H = — Yluj) XijXij + constant, where Xij = 
Su=i c l,^ c i,Ai- The model is similar to the 2D SU(N) t-J model of Affleck and Marston pT 



with the important difference that here one fermion per site implies that each flavor of 
fermions is 1/4 filled, while in their study each flavor of fermions is close to 1/2 filled, 
the case relevant to cuprates. We consider a uniform and real mean field bond amplitude 
X — {Xij)i an d examine its stability against a generalized spin density wave (SDW) state 
with four sublattices B u , v — 1,2, 3, 4. The uniform mean field state describes a flavor liquid. 
For the spin-1/2 Heisenberg model, the uniform state was found unstable against the (tt, tt) 
SDW state However, the instabilty is related to the nesting Fermi surface at 1/2 filling. 
We expect the uniform state to become stable against the SDW state at fillings sufficiently 
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far away from 1/2, as in the present case. Minimizing the mean field energy with respect to 
the SDW order parameter m, defined so that the mean occupation number for the flavor fi 
at site i C B u , (c|„Cj jM ) = (1 — m)/4 + m8 u ^, we find that the uniform bond state is stable 
against the SDW. Thus, the mean field theory suggests the ground state of ([!]) on the square 
lattice is disordered. The flavor disordered state is found to be gapless in the fermion mean 
field theory, but from previous studies of SU{2) systems, the issue of the gaplessness needs 
to be further examined. We note that the 1/4 filled uniform mean field state is unstable 



against the commensurate flux phase of flux hc/(4e) per plaquette [p^ , p3|| . The flux phase 
is also a flavor liquid. 

We expect the model on the triangular lattice to be also disordered. The AF Heisenberg 
spin-1/2 system on the triangular lattice is believed to order in a three-sublattice 120° 



structure p6| . With orbital degeneracy, such a spin ordering is no longer favored. In Fig. 
2, we compare the estimated energies for various long range ordered states, including the 
classical SU(4:) Neel state (the same as the Ising like Neel state for both triangular and 
square lattices, see reference [|24||), the orbital polarized spin ordered state, and a valence 
bond state, with the SU (4) singlet plaquette solid state. Unlike the spin only problem, where 
the classical Neel state and the valence bond solid state are degenerate, here the plaquette 
state has much lower energy than all others. Since the plaquette state can resonate to 
further lower the energy (and becomes a flavor liquid), we speculate the ground state to be 
a resonant plaquette state with neither spin nor orbital long range order. 

It is interesting to note that a spin liquid state with two-fold orbital degeneracy may 
have already been realized in the best samples of L1N1O2 ■ The measurements of magnetic 
susceptibility, specific heat, /iSR, and NMR at low temperature show no long range ordering 
in spins, and the /zSR also shows that iVi-spins remain fluctuating even at 20 mK ||. These 
evidence strongly indicate a spin liquid ground state. In that material, a formal iVi 3+ -ion 
has spin s = 1/2 and a two-fold e g orbital degeneracy. The iVi-ions form layered triangular- 
lattices, separated by two oxygen and one Li layers, so that the interlayer coupling is weak. 
We believe the orbital degeneracy is responsible for the observed spin liquid properties, and 



Eq. (jH) may serve as a simple model to illustrate the role of the orbital degeneracy. 

We now turn to the ordered states where the SU(4) symmetry is broken. There are many 
ways to break SU(4), and here we discuss the generalized 4-sublattice Neel state shown in 
Fig. 2a. The remaining symmetry is U(l) x (7(1) x U(l). There are a total of 12 Goldstone 
modes, comprised of 4 spin-wave modes, 4 orbital-wave modes, and 4 spin-orbital-crossing 
modes. 

The n.n. correlation function can have very unusual properties in states with symmetry 
breaking. Consider for example the 4-sublattice Neel state. Let the two n.n. sites % and j 
belong to the two sublattices with (sf ) = (sj) > 0. In spin only systems, this would imply 
a ferromagnetic coupling and the correlation w(s z ) = (sfsj) > 0. In the presence of the 
orbital degeneracy, the situation can be different. Let us start with the disordered state, 
where w(s z ) < from (10). Provided the transition is continuous, this implies that w(s z ) will 
remain negative at the transition point, or close to it on the ordered side of the transition. 
To estimate the critical value of (sf), we have carried out calculations on a four-site system. 
We add a term in (1) representing alternative fields coupled to the spin and orbital, so that 
the symmetry is explicitly broken, and the four sites represent four sublattices. The ground 
state and the correlation functions are calculated numerically. We find w(s z ) < when 
< ( s i) = ( s j) < 0.18. Therefore the system can have antiparallel spin correlation while 
both spins have parallel magnetic moments. This unusual magnetic property can be tested 
by neutron scattering. 

We would like to thank I. Affleck, B. Frischmuth, N. Kawakami, B. Normand, T. M. Rice, 
H. Shiba, M. Takano, K. Ueda for many stimulating and useful discussions. FCZ wishes to 
thank Yukawa Institute, Kyoto University for its hospitality during his visit, where part of 
this work was initiated. At the completion of the work, we learned that K. Ueda et al. have 
considered model (1) with emphasis in 1-dimension ||27| . 
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FIGURES 

FIG. 1. Young tableaux for SU(4) model (1) a) in a two-site system, and b) in a four-site system. 
The dimensionality of representations is indicated for each tableau, and inside the parenthesis is 
the number of distinct representation. 

FIG. 2. Average energy per bond e of (1) for various states for the triangular lattice, a) Classical 
4-sublattice SU(4) Neel state, with each flavor shown by its orbital (dashed arrow) and spin (solid 
arrow) state, b) State with orbital "ferromagnetic" and spin AF. In this case, Eq. (1) is reduced 
to H = J2(i.j)(2si ■ Sj + 1/2). e is deduced from Ref. p6|| . c) Valence bond state. Each double-line 
represents a two-site valence bond of orbital singlet and spin triplet. Note the spin long range 
order, d) Plaquette state. Each plaquette (linked by four thick lines) represents a four-site SU (4) 
singlet. 
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